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We consider the energy dependence of the local V and CV violation in Au+Au 
and Cu+Cu collisions over a large energy range within a simple phenomenological 
model. It is expected that at LHC the chiral magnetic effect will be about 20 times 



weaker than at RHIC. At lower energy range this effect should vanish sharply at 
energy somewhere above the top SPS one. To elucidate CME background effects a 



transport model including magnetic field evolution is put forward. 



I. INTRODUCTION 

As was argued in Refs. [l-4] the QCD topological effects in intimate connection with 
axial anomaly may be observed in heavy ion collisions directly in the presence of very 
intense external electromagnetic fields due to the "Chiral Magnetic Effect" (CME) as a 
manifestation of spontaneous violation of the CV symmetry. First experimental evidence for 
the CME identified via the observed charge separation effect with respect to the reaction 
plane has been presented by the STAR Collaboration at RHIC Q| . In this paper we analyze 
the STAR data in a simple phenomenological way to estimate a possibility observing the 
CME in the larger energy range, from the LHC to FAIR/NICA energies. We also make a 
step towards a dynamical estimate of the CME background based on the nonequilibrium 
Hadron-String-Dynamics (HSD) microscopical transport approach [6] including magnetic 
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field. 



II. PHENOMENOLOGICAL ESTIMATES OF THE CME 

n 

A characteristic scale of the process is given by the saturation momentum Q s [I], so 
the transverse momentum of particles p t ~ Q s . Then the total transverse energy per unit 
rapidity at mid-rapidity deposited at the formation of hot matter is expressed through the 
overlapping surface of two colliding nuclei in the transverse plane 5* 

^ ~ e-V = e.Az.S = Q s .(Q*S)~Q s .^f^ . (1) 
ay ay 

Here the energy density and longitudinal size Az ~ At ~ 1/Q S are taken in order of 
magnitude as follows e ~ and Az ~ 1/Q S - 

For one-dimensional random walk in the topological number space the topological charge 
(winding number) generated during the time tb, when the magnetic field is present, may be 
estimated as 

n w ^^ s = ^T s .V.r B ^^^^.^Q~^, (2) 

where T$ is the sphaleron transition rate which in weak and strong coupling T$ ~ T 4 with 
different coefficients. The initial temperature To of the produced matter at time r ~ 1/Q S 
is proportional to the saturation momentum Q s , T = c Q s . At the last step of ([2]) the 
expansion time and the corresponding time dependence of the temperature are neglected. 
Since sizable sphaleron transitions occur only in the deconfined phase, the time tb in Eq. (J2]) 
is really the smallest lifetime between the strong magnetic field fs one and the lifetime of 
deconfined matter r e : 

tb = mm{f B ,r t }. (3) 

The measured electric charge particle asymmetry is associated with the averaged corre- 
lator a by the following relation [7J: 

(cos(?/> Q +ipfi- 2^ RP )) = (cos(> Q + ifjp - 2ip c )) / v 2 , c = Vi, a Vi,(3 ~ a a a p , (4) 

where ^ rp is the azimuthal angle of the reaction plane defined by the beam axis and the 
line joining the centers of colliding nuclei. Averaging in (jlj) is carried out over the whole 
event ensemble. The second equality in (jlj) corresponds to azimuthal measurements with 
respect to particle of type c extracted from three-body correlation analysis Q, V\ and v 2 



are the directed and elliptic flow parameters, respectively. According to Ref. [l] an average 
correlator a = ^/a a ap is related to the topological charge, 

Tl w 

s T B l T B , ^ \— 1/16 



a ~ 



dN h 



"hadrons/ ' dy ^ dN hSudmns / dy 
where absorption and rescattering in dense matter responsible are neglected for the same 



(5) 



and opposite charge correlations. In the last equality we assumed that Q 2 S ~ s 



dN^Amns/ dy [8]. Our susequent consideration is based on Eq. (0). 

Thus, the direct energy dependence of average correlator is comparatively weak. Results 
of dynamical heavy-ion calculations of the magnetic field at the central point of the transverse 
overlapping region of colliding nuclei and energy density of created particles are presented 
in Figs. [H and El respectively. Here for a field estimate we follow Ref. [9( basing on the 
UrQMD model [10] and applying the electromagnetic Lienard-Wiechert potential with the 
retardation condition for the magnetic field. As is seen, at the impact parameter b = 10 fm 
the maximal strength of the dominant magnetic field component B y (being perpendicular 
to the reaction plane) is decreased in Au+Au collisions by the factor of about 10, when 
one proceeds from ^/snn =200 GeV to E[ a b =11 GeV, while for the created particle energy 
density e in the central box this factor is 250, i.e. noticeably higher. 
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FIG. 1: The time evolution of the magnetic field strength eB y at the central region in Au+Au 



collisions with the impact parameter b = 10 fm for different 



carried out. within the I'rOMD model [l(| (lor a detail see 9|). 



bombarding energies. Calculations are 



To use Eq. (jHJ) we need to identify the impact parameter, saturation momentum and 
multiplicity at a specific centrality. These can be found in Ref. js| where the Glauber calcu- 
lations were done. As a reference point we choose b =10 fm in our subsequent consideration. 
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FIG. 2: The time evolution of the energy density e of created particles in the Lorentz-contracted 
box with the 2 fm side at the central point of overlapping region. The impact parameter b = 10 
fm. 



The measured value of (cos(ijj a — for the same charge particles from Au+Au 

(■\/snn = 200 GeV) collisions at the impact parameter b =10 fm (40-50% centrality interval) 
is —(0.312 ± 0.027) • 10~ 3 j^j]. Appropriate number for ^snn = 62 GeV seems to be a little 
bit larger but for Cu+Cu collisions the effect is definitely stronger jsj. Thus, ignoring any 
final state interactions with medium, assuming a a = ap = a and neglecting the directed 



flow vi a = vib = we get from Eq. (TjJ a 2 = 0.31 ■ 10 3 for the maximal RHIC energy. 



exp 



Using numbers for the a/s/vaF =200 GeV reference case, from Eq. dSJ) we may quantify the 
CV violation effect by the correlator 

a 2 = K Au (v^)- 1/8 -r B . (6) 

The normalization constant Ka u can be tunned at the reference energy tJsnn =200 GeV 
from the inverse relation and experimental value a exp at this energy for b =10 fm 

_ aL p -(200) 1 / 8 
Au ~ r B (200) • {i) 

The lifetime tb may be denned as the time during which the magnetic field is above the 
critical value needed to support a fermion Landau level on the domain wall eB cr u = 2ir/Sd, 
where is the domain wall area. Since the size of the domain wall is not reliably known, 
it is hard to pin down the number, but it should be of the order of m 2 . Honestly, we have 
to treat it as a free parameter. 
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TABLE I: Estimated parameters for the CV violation effect in Au+Au collisions at centrality 
(40-50)% with the critical field eB crit = 0.2 m£. 
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1/16 
S NN 
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Indeed the size of the topological defect (say, a sphaleron) in the region between T c and 
2T C is very uncertain. At weak coupling, the size is determined by the magnetic screening 
mass and it is ~ l/(a s T). If one plugs a s ~0.5 and T = 200 MeV, the size is of about 2 fm 
and then the threshold field is very small eB y ~ (a s T) 2 ~ 0.2 m 2 . 

On the other hand, we know that between T c and 2T C the magnetic screening mass which 
determines the size of the sphaleron is not small as expected from the perturbative theory, 
a> s T, but from the lattice it is numerically large till about 5T C . This would increase the 
threshold to 20 m 2 , however the relation between magnetic mass and the sphaleron size is 
valid only as long as the coupling is weak. All we can say it is perhaps in between (0.2—20) 
m 2 . Eventually lattice QCD calculations may clear this up. 

The upper bound on the magnetic strength eB crit = 20 m 2 results in t# = even for 
the RHIC energy and therefore in this case the CME should not be observable at all in 
this energy range. The time evolution of the magnetic field and energy density, e, of newly 
created hadrons are presented in Figs. [1] and [2j The extracted values of tb defined by the 
constraints eB y > 0.2 m 2 and t € (e > 1 GeV/fm 3 ) are summed in Tabl. [H For the reference 
energy and the minimal magnetic field constraint we have K^ u = 2.52 • 10~ 3 . If lifetimes 
are known for all energies one can estimate the CV violation effect through the a 2 excitation 
function. 

From the first glimpse as follows from Tabl. [U in the case of eB cr u = 0.2 m 2 the interaction 
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time t b is defined solely by evolution of the magnetic field since f B < t £ whereas r £ rs 2 
fm independent of ^snn- The expected CME for Au+Au at b = 10 fm (see the last 
column in Tabl. H]) monotonously increases when ^/snn goes down but then sharply vanishes 
exhibiting a shallow maximum in the range between near the top SPS and NICA energies. 
The position of CME maximum and its magnitude depend on the cut level which just defines 
¥g. The decrease of the eB y bound till 0.02 m 2 shifts the maximum toward lower energy 
^snn and enhances its magnitude. In an opposite limit when results are extrapolated to 
the LHC energy, the CME falls down by a factor of about 20 with respect to the RHIC 
energy. This result is quite understandable. The CME is mainly defined by the relaxation 
time of the magnetic field which is concentrated in the Lorentz-contracted nuclear region 
~ 2R/j. Therefore, the CME is inversely proportional to the colliding energy, ~ 1/\/ Saw > 
and proceeding from the RHIC to LHC energy we roughly get the suppression factor about 
4.5/0.2^ 22. 

There is one worrying point here. Proceeding from ^snn =200 to 62 GeV the predicted 
value of a 2 for 6=10 fm increases in three times though not more 20% growth has been 
observed in these collisions in the recent experiment jsj . This essential disagreement cannot 
be removed by a simple variation of eB crit . One may try to explain this correlator overesti- 
mation at ^/snn =62 GeV by an irrelevant choice of the energy dependence of multiplicity 
in Eq. ([6]). For the correlator ratio at these two energies we have 

0^200) = / JB V» = e „ 

a 2 (62) r s (62) V2007 K ' K ' 

where we use lifetime values from Tabl. [J and experimental values for correlators Q , j3 = 
1/8. As follows from Eq. (jSJ), to explain the experiment the exponent should be negative, 
< 0. Therefore, the fast growth of t b with the energy decrease cannot be compensated 
by uncertainty in the energy dependence of the correlator a. 

Uncertainty in the choice of the impact parameter does not help us to solve this issue. 
It turned out that one fails to fit this ratio by the variation of only eB crit . Here we should 
remember that not only the strong magnetic field but also high density of soft equilibrium 
quark-gluon matter are needed. Equilibration requires some finite initial time ti iE which we 
associate with the moment when a maximum in the e(t) is achieved (see Figj2j) . This makes 
f shorter and in combination with eB crit variation, t b = f B — tj >e , allows us to satisfy the 
condition (jSJ). Using the value of r B (200) obtained in this analysis one can recalculate the 
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coefficient in Eq. ([6]), K^ u = 6.05 • 10 -3 , and therefore find the correlator a at any energy. 
In principle, similar analysis may be repeated for other impact parameters to consider the 
6-dependence of the CME. As was shown in Refs. ^] the CME roughly is linear in b/R. 
Taking this as a hypothesis we evaluate the centrality dependence of the CME fitting this 
line to points b = 10 fm (or centrality (40 — 50)%) to be estimated in our model and 6 = 
where the CME is zero. The results are presented in Fig. [3] for Au+Au collisions at three 
energies. 
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FIG. 3: Centrality dependence of the CME. Experimental points for Au+Au and Cu+Cu collisions 
are from {jj. The dotted line is our prediction for Au+Au collisions at the LHC energy. 



As it is seen the calculated lines quite reasonably reproduce the measured points of 
azimuthal asymmetry of charge separation for Au+Au collisions at a/J/vtv =200 and 62 
GeV. The chosen value of eB crit = 0.7 results in absence of the CME at the top SPS 
energy because the critical magnetic field practically coincides with the maximal field at this 
bombarding energy (see Fig. [TJ. The CME at the LHC energy is expected to be less by a 
factor of about 20 as compared to that at the RHIC energy. Note that at the LHC energy 
we applied a simplified semi-analytical model [9] for magnetic field creation and assumed 
tj j£ = 0. Thus, we consider this LHC estimate as an upper limit for the CME. 

Similar analysis can be repeated for Cu+Cu collisions basing on available RHIC mea- 
surements at two collision energies. Here one remark is in order. An enhancement of the 
CME in Cu+Cu collisions with respect to Au+Au ones was seen experimentally at the same 
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centrality j^J but not at the same impact parameter. As follows from the Glauber calcula- 
tions, the impact parameter b=10 fm for gold reactions corresponds to centrality (40-50)% 
while the same centrality for copper collisions matches b=4.2 fm. The time distributions 
of the magnetic field and energy density for Cu+Cu collisions look very similar to that for 
Au+Au ones but lifetimes, both re and r e , are shorter in the Cu+Cu case. For the extracted 
lifetimes and other characteristics at eB crit = 0.2m 2 [Kcu — 6.34 • 10~ 3 ) we meet again the 
same problem: one should compensate a too strong energy dependence of the model cor- 
relators by the proper definition of lifetimes. Defining the lifetime in the same manner as 
for Au+Au collisions the lifetime ratio tb(62)/tb(200) is turned out to be very close to 
experimental one at eB cr i t =0.3 m 2 . In this case Kc u = H-9 • 10~ 3 . In the linear approx- 
imation with the reference point at b =4.2 fm, one may draw the centrality dependence of 
the CME for Cu+Cu collisions shown also in Figj3] which is in a reasonable agreement with 
the experiment. Note that eB crit =0.3 m\ which is slightly above the maximal magnetic 
field at y/SNN =62 GeV implies that the CME for Cu+Cu collisions will not be observable 
even at the top SPS energy. 

From dimensionality arguments the system-size dependence of the chiral magnetic effect 
(at the same all other conditions) would be expected to be defined by the surface S = Sa(&) 
of an "almond" -like transverse area of overlapping nuclei since both the high magnetic field 
and deconfined matter are needed for this effect. The magnetic field was evaluated in the the 
center of the chapping region but as was shown in Re, fl the stndied eB v _nt is 
quite homogeneous along x of this "almond" . Using for "almond" area a rough estimate as 
two overlapping discs of radius R = r A 1//3 , namely 5* = S^ft) = t^\J R 2 — (b/2) 2 (R — 6/2), 
we have Scu(b = A.2)/S^ u (b = 10) ~ 1.65 which seems to be consistent with experimental 
ratio of the CME at ^/snn for these two points. However, this result was obtained for 
different eB crit and non-zero initial time t itE , and this success cannot be repeated for Cu+Cu 
(62 GeV) collisions. Therefore, the Cu enhancement effect is not only a geometric one. 

III. TOWARDS A KINETIC APPROACH TO THE CME BACKGROUND 

The discussed CME signal, the electric charge asymmetry with respect to the reaction 
plane, may originate not only from the spontaneous local CP violation but also be simulated 
by other possible effects. In this respect it is important to consider the CME background. 
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We shall do that considering a full evolution of nucleus-nucleus collisions in terms of the HSD 
transport model [&] but including formation of electromagnetic field as well as its evolution 
and impact on particle propagation. 

Generalized on-shell transport equations for strongly interacting particles in the presence 
of magnetic fields can be written as 



+ (Vjr U) V r - - (V T , U + gv x B) V p } f(r,p, t) = I coll (f, A, ...f N ) 

which are supplemented by the wave equation for the magnetic field whose solution in 
the semi-classical approximation for point-like moving charges is reduced to the retarded 
Lienard-Wiechert potential used above The term U ~ Re(Tj ret )/2p is the hadronic 
mean-field. 




FIG. 4: Distribution of the magnetic field strength eB y in the y = plane at t =0.05 and 0.2 (in 
the middle) fm/c for Au+Au collisions at ^/snn =200 and b =10.2 fm. The collision geometry is 
projected on x — z plane by points corresponding to a particular spectator position. Curves (and 
their projections) are levels of the constant eB y . 

One should note that the off-shell HSD transport approach is based not on the Boltzmann- 
like transport equation (Q but rather on the off-shell Kadanoff-Baym equations having 
similar general structure. The set of equations was solved in a quasiparticle approximation 
by using the Monte-Carlo parallel ensemble method. To find the magnetic field a space 
grid was used. In a lattice point of this grid the retarded vector potential is evaluated. 
The magnetic field is calculated by its numerical differentiation. The field inside a cell is 
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approximated by that at the nearest grid point. To avoid singularities and self- interaction 
effects, particles within a given cell are excluded from procedure of the field calculation. 

An evolution snapshot of the magnetic field B y (x,y = 0, z, t) (in units of m 2 ) formed 
in Au+Au (200 GeV) peripheral (b =10.2 fm) collisions are given in FigJH for two time 
moments t =0.05 and 0.20 fm/c. The collisional geometry is presented by a set of points 
every of which corresponds to a spectator nucleon. The whole field is not homogeneous 
exhibiting a wide maximum over the transverse size of overlapping (participant) matter and 
strong contraction in longitudinal direction. Opposite rotation of the magnetic field along 
direction of two colliding nuclei results in corresponding two minima from outer sides of 
spectator matter remnants. At expansion these remnants are moving away from each other. 
The position of a maximum in the magnetic field strongly correlates with that in the energy 
density of created particles. Large local values of B y and e reached in these Au+Au collisions 
provide necessary conditions for observation of signals of a possible parity violation. 



IV. DISCUSSION AND CONCLUSIONS 



Summarizing one should note that for heavy-ion collisions at ^/snn <i 11 GeV the mag- 
netic field and energy density of deconfined matter reach very high values which seem to be 
high enough for manifestation of the Chiral Magnetic Effect. However, these are only neces- 
sary conditions. To estimate a possible CME a particular model is needed. For the average 
correlator our qualitative prediction a 2 ~ has a rather small exponent but nevertheless 
it is too strong to describe observable energy behavior of the CME. This model energy de- 
pendence can be reconcile with experiment 5J by a detailed treatment of the lifetime taking 
into account both the time of being in a strong magnetic field and time evolution of the 
energy density in the QGP phase. For the chosen parameters we are able to describe data 
for Au+Au collisions on electric charge separation at two available energies. We predict 
that the effect will be much smaller at the LHC energy and will sharply disappear near the 
top energy of SPS. Coming experiments at the Large Hadron Collider and that the planned 
Beam Energy Scan program at RHIC [12] are of great interest since they will allow one to 
test the CME scenario and to infer the critical magnetic field eB cr i t . 

The experimentally observed CME enhancement for Cu+Cu collisions is related with the 
selection of different impact parameters for the same centrality. However, it is not reduced 
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to a purely geometrical effect. 

The problem of parity violation in strong interactions and the related CME are actively 
debated now. It is of great interest that the electric charge asymmetry with respect to the 
reaction plane may originate not only from the spontaneous local CP violation but also be 
simulated by other possible effects. First step in study of dynamical study of the CME 
background has been made in Sec lIIII It is important that the developed kinetic approach 
in principle allows one to simulate the Chiral Magnetic effect itself. This work is in progress. 
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